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A molecular description for magic-number configurations of interacting electrons in a quan- 
tum dot in high magnetic fields developed by one of the authors has been elaborated for four, 
five and six electron dots. For four electrons, the magic spin-singlet states are found to alter- 
nate between two different resonating valence bond (RVB)-like states. For the five-electron 
spin-polarized case, the molecular description is shown to work for the known phenomenon 
of magic-number sequences that correspond to both the A^-fold symmetric ring configuration 
and a {N — l)-fold symmetric one with a center electron. A six-electron dot is shown here 
to have an additional feature in which inclusion of quantum mechanical mixing between 
classical configurations, which are deformed and degenerate, restores the A-fold symmetry 
and reproduces the ground-state energy accurately. 
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I. INTRODUCTION 

Quantum dots are two-dimensional (2D) electron 
systems confined in a small few hundred nanome- 
ter) region and have attracted much attention as ar- 
tificial atoms |l] . A remarkable feature of these 
'atoms' is their controllability. A dramatic feature 
in fact arises when quantum dots are placed in high 
magnetic fields, where specific total-angular momen- 
tum (L) states are stabilized with downward cusps 
in the energy against L. The cusps, originally found 
from numerical calculations Q| g, are called 'magic 
numbers' in the angular momenta. When the mag- 
nitude of magnetic field is varied, the true ground 
state changes from one magic-number state to an- 
other. Some of the cusps occur at the same frac- 
tional filling factors as the fractional quantum Hall 
effect in bulk systems. 

In a different light, quantum dots may be regarded 
as a 2D 'molecule' that consists only of electrons, 
as contrasted with the ordinary atoms where elec- 
trons are bound to nuclei. The quantum states in 
strong magnetic fields can in fact be treated with 
a molecular description, which has been developed 
by one of the present authors [D . In this approach 
we start from the classical equilibrium configura- 
tion and make an expansion of the Hamiltonian to 
derive a quantum mechanical rotational- vibrational 
Hamiltonian to obtain anti-symmetrized wavefunc- 
tions. In a previous work this approach has been 
tested for ground and excited states of three elec- 
trons, and was found to give energies as accurate as 
within 0.1% of the exact-diagonalization result for 
L ~ 30 in a parabolic confining potential. 

Here we extend the theory to four, five, and six 
electrons in search of features for larger number of 
electrons. All our results are for a GaAs dot with a 



typical experimental situation of 4meV confinement 
energy in a magnetic field of 20T. 



II. MAGIC NUMBERS 

The magic numbers are known to depend on the 
number of electrons, N, and total spin, S. For fully 
spin-polarized electrons, magic L's are 3,6,9,12,... 
for iV = 3 and L = 6,10,14,18,... for A^ = 4. 
Namely, the magic numbers occur with an interval of 
N, or more precisely L = N{N — l)/2 + Nx integer. 
Another sequence of magic numbers, having an in- 
terval A^-1 with L = A^(A^-l)/2 + (A^-l) X integer, 
[H @ appears for six spin-polarized electrons. 
For five spin-polarized electrons, the main sequence 
has an interval A'' and the ground state for any the 
magnetic field belongs to the sequence. If one looks 
more closely at the L- dependence, however, there 
are weak cusps with an interval A^ — 1, and we can 
call these states excited magic-number states. It is 
known that an A^-fold ring has the magic number 
states with interval A^, while an A^ — 1-fold ring with 
a particle at the center has an interval A^ — 1. If we 
look at the quantum mechanical density correlation 
function in Fig. |l|, where we fix the position of one 
electron to plot the density of others, we can indeed 
identify both A^- and A^ — 1-fold rings. 




























1 



FIG. 1. Quantum mechanical density correlation with 
one electron fixed at the position indicated by a dot 
(upper panels) as compared with classical equilibrium 
point-charge configurations (lower panels) for typical 
magic-number states. For N = 5, 6, A^-fold and 
(A'^ — l)-fold symmetric states occur and give rise to two 
sequences of magic numbers with intervals A*' and TV — 1, 
respectively. 



III. ECKART-FRAME THEORY 

The molecular description provides not only a way 
of understanding the magic numbers, but also a 
quantitatively accurate description of the ground- 
state energy in the strong field regime [D . The start- 
ing point is the remarkable correspondence between 
the quantum mechanical configurations (upper pan- 
els of Fig.l) and the classical equilibrium configu- 
rations (lower panels) of point charges interacting 
with the Coulomb repulsion. The physical picture 
of the system in the high field regime is then that of 
an 'electron molecule' which rotates with total an- 
gular momentum L and vibrates around its classi- 
cal equilibrium configuration. We can first harmon- 
ically expand the Hamiltonian around the classical 
configuration. The vibration is treated in a similar 
way to that of chemical molecules, and this is facili- 
tated by the use of a molecular-fixed reference frame 
called the Eckart frame in which the rotational and 
vibrational motions decouple to first order ||]. We 
then diagonalizc the Hessian matrix for the harmon- 
ically expanded Hamiltonian and obtain the effective 
potential. Using the effective potential and taking 
account of Coriolis coupling, we find the rotational- 
vibrational states. The final important ingredient is 
the anti-symmetrization of the wavefunction. 

When the classical equilibrium configuration has 
Cm symmetry, a group-theoretical consideration dic- 
tates that antisymmetrization enforces a wavefunc- 
tion to have a factor 
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where ni,ky{i) are the quantum number and phase 
parameter corresponding to the ith normal mode, 
respectively, while kg is the phase parameter for the 
spin wavefunction. It follows that the wavefunction 
vanishes unless 
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This criterion reproduces the correlation between 
the N- and {N — l)-fold rings and the interval of 
magic numbers. 

After quantum numbers that satisfy this condition 
have been found, the total energy can be written in 
the form 



E = 



Eel + E, 



zp 
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where E^^^ is the energy of center-of-mass motion, 
Ed is that of the classical configuration, E^p is the 
quantum zero-point energy and Ex{ni, ■ • ■ , J^aw-a) 
the excitation energy of the vibrational state. 



IV. FOUR ELECTRONS 

We have calculated energies and density correla- 
tion functions for ground states against total angular 
momentum for all the spin configurations using both 
exact diagonalization and Eckart-frame theory. The 
energies for spin-singlet {S — 0) states are plotted in 
figure 1^, which shows that the Eckart-frame theory 
produces accurate ground state energies. We have 
also checked that ground-state energies obtained by 
Eckart-frame theory for all the other spin configura- 
tions agree very well with those of the exact diago- 
nalization. One can readily see from equation |^ that 
the interval between magic numbers is affected by 
the total spin S. For the spin-polarized case, 5 = 2 
here, the magic numbers are L — 6,10,14,18,... 
and these are explained by equation ^ with fc^ = 0. 
In the case of 5* = states magic numbers are 
L 6, 8, 10, 12, 14, 16, 18, .... The spin-resolved cor- 
relation functions are displayed in fig. H, which show 
that there are two different sequences; one has L = 
10, 14, ... while the other has L ~ 12, 16, .... In the 
Eckart-frame theory, these magic numbers are char- 
acterized by different kg. The former corresponds to 
fcs = and the latter to fc^ = 2. Except for a nor- 
malization, the spin wavefunctions are written as 

fcs-0: 1 TTii> + UTTi> + UiTT> + l Tiit) 

-2[| TiTi) + UTiT)], (4) 
fc. = 2 : 1 TTii) - I iTTi) + I iiTT) - I Tiit). (5) 

We can rewrite these in terms of spin-singlet pairs 
as 



(6) 
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where 
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(I Ti> - I it))- These wavefunc- 



tions are the same as the RVB (resonating valence 
bond) states considered for the Hubbard model 
The configuration of the spin resolved density cor- 
relation in Fig. 1^ is indeed explained from equa- 
tions |J|. 
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FIG. 2. (a) Lowest energies against the total angular 
momentum L for A'' = 4, 5 = 0. 

(b) The same for N — 5, S — 5/2. Solid (open) arrows 
indicate the magic numbers with interval TV (TV — 1). In 
each figure, the solid line with filled circles represents the 
result of Eckart-frame theory while the dotted line with 
open circles gives the exact diagonalization result. 
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FIG. 3. Spin-resolved density correlation functions 
for TV = 4, S = magic-number states. We fix an 
up-spin at the position indicated with a filled circle to 
contour-plot the density of parallel spins (upper panels) 
and anti-parallel spins (lower panels). 



V. FIVE ELECTRONS 

In the case of five electrons, there are two different 
classical equilibrium configurations: a pentagon and 
a square with one particle at the center. These con- 
figurations are shown in Fig. ^ We have calculated 
lowest-energy states against L in the Eckart-frame 
theory for both of the configurations, and have plot- 
ted the lowest of these two energies in Fig. S. The 



difference from the exact- diagonalization result is 
less than 0.2 % for L ~ 40. 

In Fig. ^ we have indicated the magic numbers 
with the interval TV (solid-down-arrows) and the in- 
terval TV — 1 (open-up-arrows). If we look at in Fig.l 
the density correlation for the states belonging to 
the TV-sequence {L = 30) or the (TV — l)-sequence 
(L = 22), the shape of the molecule does indeed 
possess the TV-fold or (TV — l)-fold symmetries, re- 
spectively. 



VI. SIX ELECTRONS 

For the six electron system, we have again two 
kinds of classical equilibrium configurations: a six- 
fold ring and a pentagon with a center particle as 
shown in Fig.l. If we more closely look at the 
six-fold ring, the classical configuration is distorted 
from the regular hexagon and there are doubly- 
degenerate three-fold equilibria, which would naively 
imply magic numbers with an interval of three. 

The exact-diagonalization result for the magic 
numbers for six spin-polarized electrons, however, 
consists of a series having period of six (L = 
21,27,33) along with another having period of five 
(L = 25,30,35,40,45,...). The period of six may 
at first seem inconsistent with the three-fold clas- 
sical configurations. However, this discrepancy can 
be resolved by noting that the six-fold symmetry 
is restored when a quantum mechanical mixture of 
the rotational-vibrational states centered at the two 
three-fold configurations is used to construct the ap- 
proximate electron state. In Fig. ||, the classical en- 
ergy of a rotating configuration (L = 33, no vibra- 
tion) is plotted against the distortion (see Fig. ^ B) 
that couples the three-fold rings. The saddle point 
is the regular hexagon while the two minima are the 
rings distorted in staggered directions. The height 
of the barrier (about O.OOOlmeV) is much smaller 
than the typical lowest quantum vibrational energy 
(about 1.0 meV), so that they are almost completely 
mixed to result in a restored six-fold symmetry. 

To calculate the energy of the hexagonal states we 
approximate the vibrational potential energy. A six- 
electron molecule has nine in-plane normal modes 
and the lowest has the same displacement pattern as 
distortion shown in Fig. IJ. The vibrational potential 
depends on L and is almost parabolic for L > 15, 
but for the normal mode shown in Fig. ^ there is 
some small extra structure which originates from the 
double minimum in the energy without vibration. 

There are two relevant normal modes that are 
Coriolis-coupled, and, because the curvature of the 
vibrational potential obtained by diagonalizing the 
Hessian matrix around the saddle point is small, one 
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of the solutions of the coupled equation of vibra- 
tional motion becomes unstable. However, this way 
(i.e., use of the Hessian matrix) of evaluating the 
potential is only valid in the immediate vicinity of 
the saddle point, while the wavefunction has a much 
larger extension, over which we should approximate 
the vibrational potential (inset in Fig. ^ a) by a 
parabola. We can then fit the curvature of the vi- 
brational potential in the region of about 20% of the 
mean radius, which is a few percent larger than the 
curvature obtained by the Hessian matrix, and in- 
deed makes the vibrational motion stable. This is 
how we find the ground state energy in the six-fold 
case. 

For the five-fold configuration with a center elec- 
tron, we have obtained the ground-state energy in 
the same way as for N — 3, 4. We compare the 
ground- state energies of six-fold and five-fold con- 
figurations at each value of L. In Fig. ^ we have 
plotted the lowest of the two possible energies. The 
result agrees very well with the exact diagonalization 
and the difference between the ground-state energies 
is 0.4% or less around L = 60. 
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VII. CONCLUSION AND DISCUSSIONS 

In this paper, we have confirmed the accuracy of 
the electron-molecule description of interacting elec- 
trons in a quantum dot in high magnetic fields for 
four, five and six electron systems. For four spin un- 
polarized electrons, there are two sequences of magic 
numbers that correspond to two RVB-like spin con- 
figurations. For five electrons, accurate ground state 
energies are obtained by considering both iV-fold 
and {N — l)-fold molecules. For six electrons, we 
have shown that the quantum mixing between states 
localized on classically distorted configurations leads 
to the symmetry expected from the exact diagonal- 
ization results. 
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FIG. 4. (a) The classical potential-energy surface for 
N — 6, S = 3, L = 33, plotted against the distortion 
in the lowest normal mode. Double minima correspond 
to doubly-degenerate, deformed classical equilibria indi- 
cated in insets A and C, while the saddle point B is 
a regular hexagon. In the inset of (a), the vibrational 
potential energy is plotted against the distortional vi- 
bration. 

(b) Lowest energies against L for N = 6, S — 3. 
The solid line with filled circles is obtained in the 
Eckart-frame theory with the quantum mixing, while the 
dotted line with open circles the exact diagonalization 
result. 
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